Light storage in an optically thick atomic ensemble under conditions of 
electromagnetically induced transparency and four-wave mixing 
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We study the modification of a traditional electromagnetically induced transparency (EIT) stored 
light technique that includes both EIT and four-wave mixing (FWM) in an ensemble of hot Rb 
atoms. The standard treatment of light storage involves the coherent and reversible mapping of 
one photonic mode onto a collective spin coherence. It has been shown that unwanted, competing 
processes such as four-wave mixing are enhanced by EIT and can significantly modify the signal 
optical pulse propagation. We present theoretical and experimental evidence to indicate that while 
a Stokes field is indeed detected upon retrieval of the signal field, any information originally encoded 
in a seeded Stokes field is not independently preserved during the storage process. We present a 
simple model that describes the propagation dynamics of the fields and the impact of FWM on the 
spin wave. 

PACS numbers: 42.50.Gy, 32.70.Jz, 42.50.Md 
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I. INTRODUCTION 

The successful development of practical quantum infor- 
mation applications relies in large extent on the availabil- 
ity of high-efficiency and high-fidelity memory for quan- 
tum states of photons. Recently, several promising real- 
izations of such a quantum memory were demonstrated, 
which are based on the reversible mapping of photon 
quantum states onto long-lived collective coherences in 
ensembles of atoms [ll-ll|- In the majority of these proto- 
cols, an atomic ensemble of sufficiently high optical depth 
is necessary in order to achieve high memory efficiency. 
However, in this case, some unwanted competing pro- 
cesses may interfere with quantum memory performance, 
reducing its efficiency and fidelity @-[l2- 

In this manuscript, we investigate the propagation and 
storage of weak optical signal pulses in optically thick 
hot atomic vapor based on electromagnetically induced 
transparency (EIT) d HI] . Traditionally, an EIT-based 
light storage scheme considers the simultaneous interac- 
tion of a strong control field and a weak signal fields in 
a A-type configuration, in which two ground hyperfine 
levels are linked with a common excited state, as shown 
in Fig. [Ija). In this case, the control field strongly cou- 
ples the propagation of the signal optical field with a 
collective long-lived ground state atomic spin coherence 
(spin wave) [ll, [l3] , resulting in a reduced group veloc- 
ity for signal pulses ("slow light"). Adiabatic turn-off of 
the control field maps the quantum state of the signal 
field onto the spin wave, which can be stored and later 
retrieved by restoring the control field intensity. 

Such EIT quantum memory has been realized for both 
weak classical and non-classical signal pulses (for re- 
cent reviews, see Refs. 0-@). Some recent publica- 
tions [Tsl-flTj investigated the optimal performance of 
such an EIT memory and have confirmed that high op- 
tical depth is necessary to improved storage efficiency. 
For example, 90% memory efficiency requires an opti- 




FIG. 1. (Color online) (a) The double-A system used in theo- 
retical calculations. In our case, \g) and \s) correspond to the 
*^Rb ground state hyperfine sublevels |_F, mp) = |1, 1), |2, 1), 
respectively; |e) corresponds to the excited state \F' , mpi) = 
|2, 2). Q, i}' represent the Rabi frequencies of the same con- 
trol field acting on two different transitions, while a and a' 
represent the Rabi frequencies of the signal and Stokes fields, 
respectively, (b) Sample data for light storage for signal (top) 
and Stokes (bottom) pulses at a temperature T = 70° C (opti- 
cal depth oqL = 52). Dashed lines show propagation of these 
pulses under slow light regime (constant control field). The 
black curve is a far-detuned reference pulse. 



cal depth uqL > 100 [l8|. On the other hand, an op- 
tically dense coherent atomic medium is also known to 
enhance competing undesired effects, such as resonant 
four-wave mixing (FWM) in a doublc-A configuration. 
In this FWM process, the far-dctuned interaction of the 
control field, which resonantly drives the \s) — |e) transi- 
tion [sec Fig. [TJa)], with atomic ground state coherence 
enhances the generation of an off-resonant Stokes optical 
field [a' in Fig. [2a)]. In turn, the presence of a Stokes 
field strongly affects signal pulse propagation, and the 
propagation of both signal and Stokes fields become a 
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result of the interference between regular EIT and FWM 
processes [ill, [H, [H . 

Under such conditions, the simplified treatment of an 
EIT-based quantum memory in a single-A system is in- 
complete and fails to describe light storage at high opti- 
cal depths [l^ . In this work, we explore the mechanisms 
governing propagation of signal and Stokes pulses under 
EIT-FWM conditions and develop an intuitive analytical 
treatment of this problem. In particular, we investigate 
the prospect of their light storage, i.e., a process in which 
both signal and Stokes pulses are reversibly mapped 
onto a long-lived spin coherence and thus can be faith- 
fully recreated after some storage period. Recent exper- 
iments [i^l have showed that a spontaneously-generated 
Stokes field can be detected upon retrieval of a signal field 
from a spin coherence. Based on these results, one might 
anticipate that the spin wave might function, at least to 
some extent, as a memory for both pulses. In our experi- 
ments, both input signal and Stokes fields have non-zero 
amplitudes before entering the atomic medium. We ob- 
served that under certain experimental conditions, both 
signal and Stokes pulses appear to be delayed due to the 
interaction with atoms, i.e., both signal and Stokes out- 
puts are nonzero even after the input signal and Stokes 
fields are turned off. Moreover, when the control field 
is turned off for some time and then later turned on, 
output pulses are retrieved in both channels, as shown 
in Fig. mb). However, careful experimental and theo- 
retical investigation shows that this is not a two-mode 
storage, i.e., the signal and Stokes fields cannot be stored 
independently. Instead, under EIT-FWM conditions, the 
collective ground-state spin coherence is determined by a 
particular combination of signal and Stokes fields. More- 
over, both retrieved fields are, in fact, only very weakly 
sensitive to the input Stokes field. To explain these ef- 
fects, we present an intuitive analytical picture of the 
effects of FWM on the signal and Stokes pulses and on 
the atomic spin coherence. 



II. EXPERIMENTAL ARRANGEMENTS 

A schematic of the experimental setup is shown in 
Fig. m The output from an external cavity diode laser 
(ECDL) locked to the Rubidium Di transition (A = 795 
nm) passed through a polarizing beam splitter (PBS), 
which split off a fraction of the beam for use as a reference 
frequency, while the rest of the beam was coupled into 
a single-mode optical fiber (SMF) to improve its trans- 
verse intensity distribution. All experimental data shown 
below were obtained using weak classical laser pulses. 
To ensure the best spatial overlap and mutual phase co- 
herence, all three optical fields were derived from a sin- 
gle laser beam by phase modulation. After the fiber, 
the beam passed an acousto-optical modulator (AOM), 
and the —1 diffraction beam (downshifted by 80 MHz) 
then passed through an electro-optical phase modulator 
(EOM) operating at the frequency of the ground state hy- 
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FIG. 2. (Color online) A schematic of the experimental ar- 
rangements. See text for abbreviations. 



perfine splitting of ^'^Rh [AM/(27r) = 6.835 GHz]. This 
phase modulation produced two first modulation side- 
bands at ±Ahf of nearly equal amplitudes and opposite 
phases. The zeroth order (carrier frequency) field was 
tuned to the 5^Si/2F= 2 — > 5^Pi /2F' = 2 transition and 



acted as the control field fl l2l|. The +1 modulation 
sideband, tuned near the 5^Si/2F= 1 5^Pi/2F' = 2, 
played the role of the signal field, while the —1 sideband 
acted as the Stokes field. The resulting beam was colli- 
matcd to a diameter of 1.9 mm or 2.7 mm, as specified 
below, and circularly polarized with a quarter-wave plate 
(A/4). Typical peak control and signal field powers were 
approximately 14 mW and 40 ^W, respectively. For some 
experiments, we attenuated the amplitude of the Stokes 
field as described in Refs. [121. [T5|. 

A cylindrical Pyrex cell, of length 75 mm and diam- 
eter 22 mm, containing isotropically enriched ®''Rb and 
30 Torr Ne buffer gas, so that the pressure-broadened op- 
tical transition linewidth was 27/(27r) = 290 MHz [221, 
was mounted inside tri-layer /i-metal magnetic shielding, 
in order to reduce the effects of stray magnetic fields. The 
temperature of the cell, and correspondingly the concen- 
tration of Rb in the vapor phase, was controlled using 
a bifilar resistive heater wound around the inner-most 
shield layer. Experimental temperatures ranged from 
50° C to 80° C; which corresponded to changes in Rb den- 
sity from 1.1 X 10^^ cm~'^ to 1.2 x 10^^ cm""^, and to a 
range of optical depths aoL between 10 and 110. After 
the cell, the light beam was linearly polarized with a A/4 
plate, recombined with the unshifted reference beam, and 
sent via a multi-mode optical fiber to a fast photodetec- 
tor (PD). The beat note signals between each of the -1-1 
and — 1 modulation sidebands and the reference field were 
measured using a microwave spectrum analyzer. Because 
of the 80 MHz frequency shift introduced by the AOM, 
the different beat note frequencies of each sideband with 
the reference field allowed for independent measurement 
of their amplitudes. 

Simultaneously-programmed modulation strengths of 
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FIG. 3. (Color online) Dependence of retrieved signal (blue) 
and Stokes (red) pulse energies as a function of storage time 
at r = 70° C (aoL = 52). Here, Ts = 300 fis. We normalized 
the memory efficiencies so that the zero-storage-time memory 
efficiency is unity. Experimentally, we were not operating 
under optimal storage conditions [l5l ]. and the zero-storage- 
time memory efficiency was f» 40% for the signal pulse and 
~ 5% for the Stokes pulse. 



the AOM and EOM allowed for independent control over 
the temporal envelopes of the control field and the sig- 
nal/Stokes fields. In this set of measurements, we used 
constant levels of the control field for both writing and re- 
trieval stages and truncated Gaussian waveforms for the 
signal and Stokes fields. Each slow and stored light mea- 
surement was preceeded by a 400 ^s pulse of the control 
field at maximum intensity to ensure optical pumping of 
the atoms in the interaction region into state \g). 

Similar to Refs. [l^, [HI, we extracted the spin wave 
decoherence time Ts by measuring the reduction of the 
retrieved pulse energies in both the signal and Stokes 
channels as a function of storage time and fitting to an 
exponential decay, e~*/'^=. Fig. [3] presents a sample mea- 
surement of the decay rate at T = 70°C {aoL = 52), 
for which we measured = 300 /is. To aid in compar- 
ison, we normalized the signal and Stokes retrieval en- 
ergies to their respective zero-storage-time values, which 
in this case was ~ 40% for signal and w 5% for Stokes. 
Both the exponential trend and the correspondence be- 
tween the data obtained with the signal and Stokes chan- 
nels are representative of all experimental temperatures, 
although we found that the spin-wave decoherence rate 
does have an optical depth dependence. 



III. THEORETICAL DESCRIPTION 

A. Coupled propagation of signal and Stokes fields 
in a double-A system 

In this section, we review the relevant theoretical de- 
scription of the EIT-FWM process in a double-A interac- 



tion configuration. A single A link consisting of a strong 
control field (Rabi frequency il) and a weak signal field 
(Rabi frequency a) is usually considered in the context 
of light storage under EIT conditions. In such a single A 
system under EIT conditions (i^ — fc ~ Wgg, where i' and 
Vc are signal and control field frequencies, respectively), 
the control field enables strong coupling between the sig- 
nal field and the long-lived atomic spin coherence that 
is usually described using the formalism of dark-state 
polaritons Adiabatic variation of the control field 

power reversibly transfers the signal optical field into the 
spin coherence, allowing for the realization of a quantum 
memory. 

Achieving sufficiently high memory efficiency, however, 
requires operation with an atomic ensemble at high opti- 
cal depth [TB - flTj . where the effect on signal field propa- 
gation by the off-resonant interaction of the control field 
on the |g) — |e) optical transition (Rabi frequency il') be- 
comes significant and cannot be ignored. In particular, 
it results in the coherent creation of a Stokes field (fre- 
quency v' ~ — Ahf — S, Rabi frequency a') due to 
four- wave mixing (FWM). To properly account for this 
generation [1, [^, [H, Ho, [23l - l32j , it is necessary to take 
into account this off-resonant A link, which is shown in 
Fig. [Ha). 

The rotating-frame Hamiltonian describing the atomic 
response to the light fields in such a configuration is: 



H = -{S-Ss)\s){s 

- a\e){g\+n\e){s 



i5-25,)\e){e 



^hf 



s){g\+U.c. 



(1) 



Here, the Rabi frequencies of the signal and Stokes fields 
are a = Eji/h and a' = E'fx'/h, where E and E' are the 
slowly- varying envelopes of the signal and Stokes electric 
fields, correspondingly, yu and ^' are the real dipole ma- 
trix elements of the respective transitions. We define the 
optical polarization P[z,i) ~ pcg{z,t)\/N and the spin 
coherence S{z,t) = psg{z,t)^/N, where pij{z,t) is the 
appropriate slowly-varying position-dependent collective 
density matrix element and N is the number of atoms 
in the interaction volume. We use Floquet theory [3^ to 
adiabatically eliminate the off-resonant interaction via fl' 
and a'. To linear order in 1/Ahf and in a', one obtains 
an effective Rabi frequency il'a'* / Am- States |e) and 
\g) acquire small light shifts Ss = |r2'P/Ahf and —6s, 
respectively. For our Clebsch-Gordon coefficients fl5| . 
ft' = —\/3fl. In order to match the equat ions of motion 
for quantum fields, we define g^N = ^^on^cjl as the 
coupling constant between the signal field and the atomic 
ensemble. We implicitly assume that the frequencies of 
the signal and Stokes fields are approximately the same, 
and we rescale the light field envelopes by defining dimen- 
sionless light field envelopes £ — ^E and £' = j^E' . 

In the dipole approximation, assuming that at all times 
most of the atoms are in \g), and to linear order in the 
weak light fields £ and the atomic evolution and light 
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FIG. 4. (Color online) Storage and retrieval of 16 /is-long (FWHM) truncated Gaussian pulses at T = 70° C {aoL = 52), for a 
two-photon detuning of (a) S = —20 kHz, (b) S = kHz, and (c) S = +20 kHz. In all plots, the top graphs are experimental 
data and the bottom graphs are the theoretical predictions from Eqs. ((2][5]). The black curve is a far-detuned reference pulse; 
the blue (red) traces are the signal (Stokes) pulses. Dashed (solid) lines correspond to slow (stored) light experiments. In the 
theory plots, the green curve corresponds to a model consisting only of EIT, to provide contrast with the EIT-FWM model. 



propagation equations read [III [S US HI, HI] : 

{dt+cd,)£ ^ig\/NP, (2) 
{dt + cd,)£'* - -igVN-^S, (3) 

dtS^-T^S + i^lP + i-^gVNE'\ (4) 

dtP = -TP + iVLS + ig^/NE, (5) 

where we have defined Fq = — i{5 — 6s) and F = 
7 — i((5 — 25s). The polarization decay rate 7 and the 
spin decay rate 70 have been introduced. 

Equations dlM]) fully describe the propagation of the 
light fields and the dynamics of the spin wave and of the 
optical polarization during all stages of light storage. In 
the slow light regime, when the control field is constant in 
time \^{t) = ri], Eqs. (I2][5]) can be Fourier transformed, 
and Eqs. dH O can be solved analytically [ul d, HO, 
[28l . [29| . In the stored light regime, when the control 
field intensity is time-dependent, these equations can be 
solved numerically. 

Figure U displays the results of storage experiments 
with 16 /xs-long truncated Gaussian pulses at T = 70°C 
(optical depth a^L = 52) along with the correspond- 
ing theoretical predictions, which are obtained by nu- 
merically solving Eqs. (OS]) with the appropriate param- 
eters. We measured the control field power to be 4.7 
mW, and the beam diameter was 2.67 mm, which corre- 
sponded to 0/(27r) = 9.6 MHz, and induced a light shift 
of (5s = 17 kHz. The spin wave decay rate was measured 
to be approximately 300 /US, thus 7o/(27r) w 270 Hz. The 
results from the slow light experiment (dashed lines) and 
the stored light experiment (solid lines) are overlaid to 
facilitate shape comparison. 

For a small negative two-photon detuning (5 = — 20 kHz 
[Fig.m^a)], the signal field (in blue) experiences some dis- 
tortion during propagation [as evidenced by the bumps 
in the leakage portion of the pulse (when t < 0), which 



exits the cell before the control field is extinguished] , but 
the shape of the slow pulse is preserved during the stor- 
age process. Likewise, the fraction of the Stokes field 
that exits the medium at < > in the slow light experi- 
ment (dashed red) matches the retrieved Stokes field in 
the stored light experiment (solid red). There is an ex- 
cellent agreement between the experimental observations 
and the numerical model predictions. The green trace in 
the theory plot corresponds to standard EIT-based light 
storage of the signal field, where the FWM process has 
been artificially turned off. To compute this trace, the 
Stokes contribution in Eq. ^ is set to zero, and Eqs. ([2J 
m [5]) are solved numerically. We include this trace in or- 
der to showcase the effects of four- wave mixing on signal 
pulse shape and delay. 

Figure|4l[b) demonstrates the excellent correspondence 
between experiment and theory for a two-photon detun- 
ing of (5 = kHz. For this value of 5, the signal pulse 
is less distorted during propagation, but the pulse shape 
is still distinct from the barc-EIT model. Likewise, Fig. 
Sl^c) depicts the results for 5 = +20 kHz « Sg, where 
the two-photon detuning effectively cancels the light shift 
during the writing and retrieval stages. Under this con- 
dition, the signal pulse will experience the least amount 
of distortion due to FWM, since the EIT transmission 
peak is, at least for a sufficiently narrow pulse bandwidth, 
symmetric about 5 — Ss- As a result, the dispersion ex- 
perienced by the pulse is mostly linear. In all cases, the 
theoretical model matches the experimental data very 
well. 

The correspondence between slow light pulseshapes 
(dashed lines) and the shapes of the retrieved pulses 
(solid lines) illustrates an important result — that when 
the writing and retrieval control field amplitudes are con- 
stant in time, the process of switching the control field 
off and on has little effect on the signal and Stokes fields, 
apart from a delay and the spin wave decay during stor- 
age time. In this case, we can further understand the 
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effects of FWM by using the closed form solutions to the 
Fourier transformed versions of Eqs. d^HS]) [ill, [H HO, HI] . 
In the Appendix, we detail the derivation of the following 
two approximate equations, which intuitively describe 
the effects of FWM and KIT on pulse propagation for 

I 

£{L, t) w f (0, t - L/vg) + dt'e 

Jo 

£'* (L, t) « £'* (0, t) + Al[ ' dt'£'* {0,t 

Jo 



the case i5 = i5s. Although these equations make a set 
of strong assumptions, including the assumption of an 
infinitely wide EIT transmission window F^; — > oo, they 
preserve the essential physics in the limit of weak FWM. 
Defining Ar, = — fi^/Ahf, the equations are 



,t-t')t' + iAn dt'£'*{0,t-t'), (6) 

Jo 

t'){L/vg-t')-iAR[ ' dt'£{0,t~t'). (7) 



These equations clearly show how the effects of FWM 
grow with optical depth agL. The first term on the RHS 
of Eq. describes the delay that the signal field ex- 
periences during propagation in an EIT medium, where 
Vg = 2fl'^/{aoj) is the EIT group velocity [ll|. Due to 
the effects of FWM, the signal field acquires a small in- 
phase gain of order Ap(L/vg)^ ^ (aoi)^7^/A^f from 
times up to L/vg earlier. The farthest away times 
are weighted more heavily. Additionally, the signal 
field acquires an i-out-of-phase contribution of order 
lA^lL/vg ^ ao^T/Ahf from the Stokes field up to L/vg 
earlier with all times contributing equally. The Stokes 
field propagates undistorted and largely undelayed, but 
gets a small [order (aoi)^7^/A^f] in-phase gain from 
times up to L/vg earlier, with closest times weighted 
more heavily, and also an i-out-of-phase contribution 
of order aoLj / Ay^f from the signal field, with all times 
weighted equally. Notice that in both equations, in the 
regime where the first term on the RHS is large, small 
in-phase (Q;oi)^7^/A^f terms and small z-out-of-phase 
aoLj/Aiif terms contribute at the same {aoL)'^j'^/Alf 
order to the absolute value of the field (which is what 
our experiment measures). 

However, the first terms on the RHS are not always 
dominant. In particular, for t > 0, the first term on 
the RHS of Eq. ^ vanishes, in which case \£'*{L,t)\ 
is dominated by the third term with a small correction 
from the second term. This means, as we will confirm 
experimentally, that the retrieved Stokes field is largely 
determined by the input signal, and not by the input 
Stokes field. Similarly, if EIT group delay is comparable 
to the signal pulse duration, then, for t < 0, the RHS of 
Eq. © is small and \£{L,t)\ is significantly affected by 
the third term. This means, as we will confirm exper- 
imentally, that the signal pulse leakage can be strongly 
affected by the Stokes input, in contrast to the retrieved 
signal pulse, which is affected by the Stokes input only 
weakly. Eqs. (|6l[7]) also show that the perturbative treat- 
ment of the effects of FWM, employed to derive them, 
breaks down when |ARL/ug| > 1, i.e., when the optical 
depth is aoL > 2 Am/ J ~ 100. 

To test the validity of Eqs. ([HllTl), in Fig.lSJa), we com- 
pare the solutions obtained by numerically solving Eqs. 



i^SSii (sohd lines) to the predictions of Eqs. (0 (dot- 
ted lines). In the dashed traces, we include the results 
of a useful intermediate approximation, which does not 
assume infinite Fe and is described by Eqs. (|A4IIA5[|A11I 
lAlsp in the Appendix. For these plots, fi/(27r) = 10 MHz 
and aoL = 80; the pulse bandwidth was Auj = O.IFe. 
From the excellent correspondence between theoretical 
models, it is evident that the approximations made in 
deriving Eqs. (|6l [7]) arc valid. 

B. Effect of four- wave mixing on the spin wave 

While the solutions of Eqs. (t^JE]) accurately describe 
the evolution of light pulses and atomic variables un- 
der slow light and storage conditions, we have not yet 
used them to elucidate the role that the Stokes field plays 
in the creation of the spin coherence. Specifically, it is 
not yet clear whether the quantum memory description 
based on the dark state polariton principle [l3| is valid 
under EIT-FWM conditions. In what follows, we de- 
velop a more transparent description of light storage in 
a doublc-A system and show that in this case, the spin 
wave is determined by a particular combination of signal 
and Stokes fields. 

We obtain this result by adiabatically eliminating the 
optical polarization P{z,t). We set the time derivative 
to zero in Eq. ([5]) and find 

P{z,t)^^p{z,t) + ^^£{z,t). (8) 

Inserting Eq. into Eq. Q, we obtain the following 
equation for time evolution of the spin wave S{z, t): 

dtSiz, t)^- {v^ + 5(z, t) - .gViV^^(z, t). (9) 

It is easy to see that the spin wave depends only on a 
combination T{z,t) of signal and Stokes optical fields, 
defined as 

Tiz,t)^£iz,t)-i^£'*iz,t). (10) 
Ahf 
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FIG. 5. (Color online) (a) Results of a numerical investigation 
of slow light with a 6.6 /xs-long pulse (reference in black), with 
n/(27r) = 10 MHz and aoL = 80, so that the bandwidth of the 
pulse Alo — O.IFe. Blue traces are the signal field; red traces 
are the Stokes field. Solid lines are the result of numerically 
solving Eqs. (OS]). Dotted lines are the result of the infinite- 
Fe approximation in Eqs. H6I7|) . Dashed lines correspond to 
results obtained using numerical integration of expressions in 
Eqs. HA4IIA5|[XlHA15l) in the Appendix, (b) The spin waves 
created at the time corresponding to a 5 /is delay in Fig.[5ja). 
The solid black line is the result of numerically solving Eqs. 
(I2J5I). Dotted lines are the results from Eq. (|15p . Dashed lines 
are the results from Eqs. HA171IA23|) . 



Eq. (|9]) is analogous to the spin wave expression obtained 
through a similar treatment of a standard three-level 
light storage model [T^ . 



dtSiz,t) 



r, + ^)siz,t)-gVNpiz,t), (11) 



but with one modification — the single light field (sig- 
nal) is now replaced by a combined signal-Stokes field 
T. Thus, one might expect that it should be possible 
to store information about this joint mode in the spin 
coherence. However, only a small (aoi)^7^/A[^f fraction 
of the Stokes field [the second term on the RHS of Eq. 
([7])] exits the medium at i > after the input Stokes 
has been turned off, while the signal pulse is delayed in 
its entirety [the first term on the RHS of Eq. ([6])]. As a 
result, in contrast to the information encoded in the sig- 
nal field, most of the information encoded in the Stokes 
field is lost to leakage, which leaves the interaction region 



FIG. 6. (Color online) Results from numerical evaluation of 
Eqs. ((3l l9l [T4t (solid blue lines) and of the homogeneous ver- 
sion of Eq. (|14p (dashed blue lines) for a range of optical 
depths qqL, as indicated in the legends. The bandwidth of 
each input pulse was Auj = 0.05Fe ~ f^^/ (20^/00^/27). 



before the control field is shut off. 

The similarity between Eq. ^ and Eq. pT|) motivates 
a more detailed comparison of our EIT-FWM system 
with the traditional EIT configuration. The propagation 
equation for J- is easily obtained from the appropriate 
combination of Eqs. ([2J[3]): 



{dt + cd,) T = -—r - ^"-^s - 



£'*.{12) 



This equation is also similar to the signal propagation 
expression in the classic stored light model |14| . 



T 



T 



(13) 



except for the optical-depth-dependent Stokes term, 
which describes the generation of signal from Stokes 
during propagation through a sufficiently optically-thick 
medium. 

When the two-photon detuning is chosen such that the 
light shift is canceled {i.e., 5 = 5s), the propagation equa- 
tion becomes, to ©(l/Ahf), 



[dt + ccos^e{t)d^] F{z,t) sa iAR£'*(z,t) 
with the angle 9{t) given by tan^ 0{t) = ^t^- 



(14) 
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Analysis of above equations demonstrates two regimes 
for light storage under EIT-FWM conditions. At low op- 
tical depths, the contribution of the Stokes field on the 
RHS of Eqs. p^[H|) is negligible. In this case, the equa- 
tions for the joint field T and spin wave S become identi- 
cal to those for £ and S in the regular EIT configuration. 
For example, if we replace the RHS of Eq. (|14p with zero, 
it would describe the propagation of J- without distortion 
at a reduced group velocity Vg = c cos^ 6 « . How- 
ever, at low optical depths and t > (after the input 
Stokes pulse has been turned off), the contribution of 
the Stokes field into T is also negligible: it is small not 
only because of the small factor F/Aht in Eq. pUj) but 
also because £'{z, t) itself is small [since the first term on 
the RHS of Eq. ([7]), generalized to arbitrary z, vanishes 
for t > Q\. Thus, signal field propagation can be analyzed 
using a three-level single A, even though the Stokes field 
can be significantly affected by control and signal fields, 
as is evident from the dominance of the last term on the 
RHS of Eq. dZl) for t > 0. 

However, at higher optical depths, the term on the 
RHS of Eq. ([T^ becomes significant. Specifically, this 
term results in gain or loss of the signal field due to 
the Stokes field. The dashed blue lines in Fig. [S] depict 
the results of numerical calculations of the homogeneous 
form of Eq. pij) . Solid blue lines show the results of the 
numerical evaluation of the full form of Eq. (jl4p with 
Eqs. P HI). For these calculations, rj/(27r) = 8 MHz, 
j/{2tt) ~ 150 MHz, and the signal pulse was chosen so 
that its bandwidth, Auj = O.OSFe. It is evident from this 
graph that the Stokes contribution is not negligible for 
optical depths aoL > 25, when the simple slow propa- 
gation of breaks down due to the Stokes term on the 
RHS of Eq. HID. 

As shown in the Appendix, the same approximations 
that lead to Eqs. (|6l7p give the following expression for 
the spin wave S'(z, t) in the limit when Fe — >■ c»: 



gVN 



£{0,t ) + A| / dt'£{0,t-t')t 



(15) 



-^Ar / ^ dt'£'*{0,t-t') 
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In Fig. Eljb), we compare the shape of the spin wave 
that is obtained by numerically solving Eqs. ([2][5]) (solid 
lines) to the predictions of Eq. ((T5|) (dotted lines). As in 
Fig. [Sja), we also include the predictions of an interme- 
diate approximation, which does not assume an infinite 
Fe and is described in Eqs. (|A17[ IA23P in the Appendix. 
The reasonable agreement between the three curves in 
Fig.[5jb) implies that Eq. ([15]) does indeed contain the es- 
sential physics. In particular, under this approximation, 
the spin wave is proportional to the signal field only, as 
in a traditional three-level singlc-A EIT system [sec Eq. 
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FIG. 7. (Color online) An illustration of the modified storage 
description, (a) During the writing stage, the input signal 
field £ (Top) and Stokes field £' (Middle) propagate at dif- 
ferent group velocities through the atomic medium, creating 
a spin wave (Bottom), (b) During the storage stage, the con- 
trol field is turned off and no light fields are present. Some 
portion of the signal field has propagated through the cell and 
leaks out before the control field is extinguished. At the same 
time, most of the information in the Stokes field is lost in the 
leakage, since, in the regime (aoL)^/A <^ 1, the propagation 
of the Stokes field is affected by the atoms only weakly [see 
Eq. ((7|)]. The spin wave is preserved during storage, (c) Dur- 
ing retrieval, the control field is turned on, releasing the spin 
wave into both the signal and Stokes fields, which exit the 
vapor cell. 



(|A23|) in the Appendix], 



Siz,t) 



N 



£(z,t). 



(16) 



Moreover, under this approximation, £(z,t) [and hence 
S{z,t)] is mostly determined by the usual slowed-down 
version of the input signal [first term in the square 
brackets in Eq. (|15p ] with small corrections of order 
\/S.^\L/vg ^ aoL^/ Ahf [third term in the square brack- 
ets in Eq. (fTS)) ] and (aoi)^7^/^hf [second term in the 
square brackets in Eq. ()15|)]. 

Fig. [7] illustrates an intuitive way to understand stor- 
age under EIT-FWM conditions. At the beginning of 
the writing stage, shown in Fig. [71(a), the control field (in 
black) prepares the atoms and causes the input signal 
field £ (in blue, top) to propagate at a reduced group ve- 
locity. The Stokes pulse (in red, middle) enters the cell 
and is not completely extinguished inside the medium 
even after the reference pulse would have left the medium. 
A collective spin coherence is created in the atomic vapor 
cell (in green, bottom). As the pulses propagate through 
the atomic medium, as shown in column (b), they expe- 
rience mutual interference effects and may become dis- 
torted. The spin wave propagates along with the signal 
field. The contributions to the spin wave are determined 
by the joint mode J-', and we can discern between the con- 
tributions to the spin wave from the signal field (shown in 
light green) from those of the Stokes field (shown in dark 
green) [see Eq. ([15])]. In the regime of weak FWM, the 
propagation of the Stokes field is only weakly affected by 
the atoms [see Eq. ([7])], so that much of the Stokes field 
leaves the end of the vapor cell as leakage. Any informa- 
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tion contained in this leaked field is lost for the storage 
process, which commences when the control field is shut 
down. After some time [column (c)], the control field is 
turned on, and the spin wave is released into both the sig- 
nal and Stokes fields. It is important to note that, since 
the joint mode T is not a normal mode, the proportion 
of Stokes to signal is not fixed. 

In the regime of weak FWM (oiqL < 25), the joint 
mode J^{z, t) is determined mostly by the input signal 
field £. Thus, the propagation dynamics experienced by 
the signal pulse will be only slightly sensitive to the am- 
plitude of the seeded Stokes pulse [of order aoLj/Ahf, 
see the last term in Eq. and consequently the spin 
wave created will have the same weak dependence on 
the seeded Stokes field [see the last term in Eq. ([T5])]. 
Since the spin wave is only weakly dependent on the in- 
put Stokes field, it is possible to create approximately the 
same spin wave for different input combinations of sig- 
nal and Stokes fields. The retrieval from the spin wave 
into the light fields will consequently have this same weak 
dependence on the input Stokes field. 

Eqs. (|6I7[) support this conclusion. Specifically, the 
amplitude of the retrieved signal field [Eq. ([6])] is deter- 
mined primarily by the input signal field (the first term 
on the RHS) with a small (00X7/ Ahf)^ correction from 
the input Stokes field (the third term on the RHS). Sim- 
ilarly, since the first term on the RHS of Eq. ([7]) van- 
ishes for i > 0, the retrieved Stokes field is also deter- 
mined primarily by the input signal (the third term on 
the RHS) with a small {aoL'y / A^)'^ correction from the 
input Stokes (the second term on the RHS). At the same 
time, signal and Stokes outputs are more strongly af- 
fected by the Stokes input for t < (leaked pulses) than 
for < > (retrieved pulses). This statement is obvious for 
the Stokes field, since the first term on the RHS of Eq. 
(O does not vanish for i < 0. The reason this statement 
holds for the output signal is that the first term on the 
RHS of Eq. de]) is smaller for t < than for t > for a 
sufficiently large group delay, while the third term on the 
RHS of Eq. dS]) is larger for t < than for i > since 
f '(0, t ~ t') vanishes for t-t' >Q. 

In Fig. |51 we display the results of an experiment de- 
signed to test these conclusions. The top graph [Fig.|Sl[a)] 
depicts the storage of a 15 /is-long (FWHM) truncated 
Gaussian signal field at T 70° C [aoL = 52). The sohd 
trace corresponds to approximately equal amplitudes of 
input signal and Stokes optical pulses, while the dashed 
traces correspond to a reduced initial Stokes amplitude 
£'*{Q^t) = -0.555(0, i). 

Notice the difference in the leakage portion {t < 0) of 
both the signal pulse and the Stokes pulse as we go from 
solid curves to dashed curves, which exemplifies that both 
signal and Stokes outputs for t < do depend strongly 
on the amplitude of the seeded Stokes pulse, as we have 
explained theoretically above and as we have reported 
previously [l3|. At the same time, the retrieved {t > 0) 
Stokes and signal pulses arc both almost independent of 
the amplitude of the seeded Stokes field, which is consis- 



1.2-1 




-20 20 40 

time (|js) 



FIG. 8. (Color online) Storage of a 15 (FWHM) trun- 
cated Gaussian pulse at T = 70° C {aoL = 52) under differ- 
ent Stokes seeding conditions. Solid lines depict storage when 
the Stokes seed amplitude is the same as the signal amplitude. 
The dashed lines correspond to the case of a reduced input 
Stokes field. The black traces show reference (input) pulses, 
and the dashed black trace in the bottom plot illustrates the 
reduced Stokes seed amplitude. 

tent with the theoretical explanation above. We repeated 
similar measurements many times under a wide range of 
experimental conditions and found the retrieved pulses 
to be weakly affected by the seeded Stokes amplitude as 
long as the input signal field is comparable to or stronger 
than the input Stokes field. 



IV. OPTICAL DEPTH DEPENDENCE OF THE 
STOKES FIELD 

In this section, we present the results of storage exper- 
iments at increasing optical depths. 

Figure |9] depicts the evolution of the Stokes and signal 
fields under storage conditions as optical depth increases. 
The general features of these results are well-explained by 
the simple signal and Stokes expressions in Eqs. ([6][7l), as 
described below. In Figs. [9lja, a'), we show the results 
of slow light and stored light experiments using a 6 fis- 
long (FWHM) truncated Gaussian pulse at a tempera- 
ture of r = 50° C, which corresponds to an optical depth 
of aoL = 10. The signal pulse [blue trace in Fig. EJa)] 
experiences a reduction in group velocity during prop- 
agation, as seen by comparing the dashed trace (slow 
light) to the black trace, which is a far-detuned reference 
pulse. The Stokes pulse [red trace in Fig. closely 
mimics the far-detuned reference pulse, indicating that 
four-wave mixing is not a dominant process at this optical 
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FIG. 9. (Color online) Stokes behavior for increasing optical depths. For all cases, the two-photon detuning 5 = 0. The 
black trace is a reference pulse; the blue (red) trace is the signal (Stokes) pulse, (a, a') Storage and retrieval of a 6 ^s-long 
(FWHM) truncated Gaussian pulse at T = 50°, which corresponds to an optical depth of aoL = 10. Here, ^l/{2n) — 8.3 MHz. 
(b, b') T = 67°{aoL = 41), pulse duration is 6 /xs, n/(27r) = 7.1 MHz. (c, c') T = 76°{aoL = 82), pulse duration is 20 /is, 
n/{2n) = 12.7 MHz. (d, d') T = SO°{aoL = 110), pulse duration is 20 /is, n/(27r) = 7.8 MHz. 



depth. In a separate run, we investigate storage of these 
pulses by turning off the control field for 100 fis. Upon 
retrieval, the signal field [solid blue trace in Fig. [9l[a)] is 
recovered with a modest reduction in amplitude due to 
spin wave decay during the storage time, but its shape 
is preserved. Additionally, we retrieve a small pulse on 
the Stokes channel [solid red trace in Fig. IHlJa')]. Eq. © 
predicts that at low optical depths [{aoLj/ A^) <C 1] the 
retrieved signal pulse will be a delayed version of the in- 
put pulse (if one accounts for the storage time) , but with 
a slight modification due to the Stokes field [the last term 
in Eq. Likewise, the Stokes field will be mostly un- 
affected by the atoms [the first term on the RHS of Eq. 
(O], so most of it will leak out (see t < 0). However, 
a small Stokes pulse oc {aoL/Ahf) generated from the 
input Signal [the last term in Eq. ([7])] will be retrieved. 

Figures M)^- b') show the results of similar experi- 
ments at T = 67° C, corresponding to an optical depth 
of aoL = 41. Here, the signal shape is again preserved 
by the storage process. The four-wave mixing effects are 
exhibited by the Stokes gain in the leakage portion of the 
pulse [see Fig. |9l[b')], which leaves the interaction region 
before the storage stage occurs. This gain is described by 
the last two terms in Eq. (O. At this increased optical 
depth, the last term in Eq. ^ also predicts an increased 
Stokes output for t > 0. The effects of FWM are also ap- 
parent in the distortion that the signal field experiences 
during propagation. 

Figures c') depict the storage experiments at 

T = 76 {aoL = 82). We used a longer pulse (FWHM 
of 20 /is). At this optical depth, the Stokes pulse expe- 
riences more gain during propagation. Storage and re- 
trieval, however, still preserve the shapes of both the sig- 



nal and the Stokes pulse. Again, it is clear that the Stokes 
field gain predicted by Eq. ([7]) becomes more apparent at 
higher optical depths. We also see that, at aoL = 110 
[column (d)], the Stokes field amplitude is smaller than 
at aoL = 82 [column (c)]. This effect is most likely due to 
the absorption of the control field by unprepared atoms 
that enter the interaction region during the waiting time. 
We also note that, at aoL w 100, the pcrturbative ex- 
pansion used to derive Eqs. (|6l7p breaks down. 

V. CONCLUSION 

We have studied the phenomenon of stored light under 
conditions of electromagnetically induced transparency 
(EIT) and four-wave mixing (FWM) in an ensemble of 
hot Rb atoms. In particular, we have investigated the 
prospect of simultaneously storing both a signal and a 
Stokes pulse in a single atomic coherence, and have shown 
that independent storage of two modes is not possible. 
The reason is that most of the Stokes pulse leaks out 
of the medium during the writing stage, so that during 
retrieval both output fields are determined primarily by 
the input signal field and depend on the input Stokes 
field only very weakly. We presented a theoretical model 
based on a simple double-A system, which agreed very 
well with experimental observations. This model allowed 
us to derive a simple relationship between input and out- 
put fields, which explained the above mentioned impos- 
sibility of two-mode storage. Furthermore, we showed 
that a particularly convenient description of storage in 
an EIT-FWM system involves a joint signal-Stokes mode, 
whose dynamics we also studied. Quantum properties of 



10 



our system will be investigated in the future using the 
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Appendix A: Derivation of Eqs. ([6]), 1^, (|15|) . and 

<m 

In the main text, we omitted the derivations of Eqs. 
(ISl), ©, (|T51) . and (IT51) . In this Appendix, we present 
these derivations. 

Since experiments and numerics show that turning the 
control field off and back on has a negligible effect on 
the fields except for a delay and spin wave decay during 
the storage time, we solve Eqs. ([2][5]) in the main text 
assuming a constant control field. In the co-moving frame 
(dt + cdz — >■ cdz), Fourier transforming in time (i — !■ w 
and dt — iw), Eqs. I^SS^ can be written as 



where F + (F - iuj)(To - ioj) [see Eq. (2) in Rcf. 

a. 

To gain some intuition for how FWM may result in 
amplification, one can consider a simple case, in which 
the diagonal terms in the matrix Al in Eq. (jAip van- 
ish (equivalently, one could consider the case where the 
Stokes field also propagates in its own FIT medium). Ap- 
proximating further F ^ OF' , we find that 



M 



"07 



2A 



hf 



-1 

1 



(A2) 



and has eigenvectors (1, ±i) with eigenvalues zfc ,"^^ , cor- 



responding to an exponentially growing solution and an 
exponentially decaying solution. In our experiment, how- 
ever, the diagonal terms for the signal and the Stokes 
fields are very different. Moreover, the effect of FWM is 
rather small and can, in fact, be treated perturbatively, 
as we will show below. 

Wc checked numerically that the last entry in the ma- 
trix M in Eq. (|Aip does not significantly affect our re- 
sults. For example, it gives a contribution to £'[L,u) 
of order Q!oL7^/Ay., which will be negligible relative to 
other contributions of order (Q;o-L)^7^/A^f since a^L ^ 
1. We will therefore neglect the last entry in the matrix 
M in Eq. (jAip for the rest of this Appendix. 

Eq. d^T]) can then be solved to give [li|, [H, [13, [M [1^ 



£'*{Z,L0) 



cosh(^z) -I- ij sinh(^z) 



.•2Ab 



sinh(^z) 



sinh(^z) cosh(^z) - if sinh(^z) 



£'*(0,a;) 



1^7 Am, /3 = y(TW^F+4A|, a - ^(zFq + c.), and ^ = ^/3. 



where Ar 

Using the convolution theorem, we then obtain 



£{z,t)^ / dt'£{Q,t~t')fi{z,t')+ / dt'£{Q,t-t')f2{z,t')+ / dt' £'*{{),t - t')h{z,t'), 



£'*{z,t)=£'*{Q,t)+ / dt'£'*{Q,t^t')g2{z,t')+ / dt' £{Q,t - t')g:i{z,t'), 



where 



h{z,t') 


~ — f due 
2ttJ 


h{z,t') + h{z,t') 


= — / due 


h{z,t') 


= — / due 


5{t')+g2{z,t') 


= — / due 
2ti J 


g3{z,t') 





2i(TZ ^ — iojt 



cosh(^z) -I- i— sinh(^z) 
2 A 

^i_R.sinh(^z)e-*'^* , 



cosh(^;i 



i — sinh(^z) 



-iujt' 



-iujt' 



(A3) 

(A4) 
(A5) 

(A6) 

(A7) 

(AS) 

(A9) 
(AlO) 
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Here /i and /2 arc defined in such a way that /i captures 
pure EIT, while /2 describes how FWM changes the rela- 
tionship between the input signal and the output signal, 
/a describes the effect of the input Stokes field on the 
output signal. Similarly, the first term in Eq. (|A5[) de- 
scribes pure undistorted propagation of the Stokes field 
in the absence of FWM. 52 describes how FWM changes 
the relationship between the input Stokes field and the 
output Stokes field. Finally, 53 describes the effect of the 
input signal on the output Stokes field. 

To get some insight into the behavior of fi and gt , we 
consider the case S = Sg (generalization to arbitrary S is 



straightforward). Wc further take the limit 70 = 0, which 
is a reasonable approximation in our experiment, except 
during the waiting time between writing and retrieval 
(however, again one can easily generalize the derivation 
below to 7o 7^ 0). Furthermore, we expand /2 and 52 to 
second order in 1 / Ahf, and expand and 33 to first order 
in 1/Ahf; in other words, wc treat FWM pcrturbatively, 
which is a good approximation in our experiment, except 
in Figs.[9jc, d). Furthermore, we approximate |12| 2i(T — > 
~ where Vg = ^ is the EIT group velocity 

and Fe = — , is the width of the EIT transparency 

7^ aoL/2 

window. We then find 



FeC ^e — 



/2(^,0«A2 

« A^f'n[0,zK](i') 



e E4 



5{t' - z/vg), 



!^(sign[t']-^Erfir-(^/"^-*') 



2^z/L 



iARn[o,z/«<,](i'), 



52(^,0 « a; 



TIL 



^l{z/vg-t')Ti[Q,zlvg]{t'), 



gi{z,t') = -h{z,t'). 



Sign[t'] 



(All) 

(A12) 
(A13) 



(A14) 
(A15) 



Here Erf is the error function, the sign function Sign[t'] ~ 
1 for i' > and —1 otherwise, and the box function 
H[a;, j/](t) = \iox X <t <y and otherwise. The second 
approximation in Eqs. (jAlHIAlS]) is done in the limit 
Fe — >■ 00 (the case of an infinitely wide EIT window). 
Using the Fe — )■ 00 expressions, we arrive at Eqs. ([HIE]). 

In Fig. [TOla-d'). we plot functions /, and for j = 
1,2,3 and the two approximate forms described above. 
Red curves depict the results of numerical integration of 
Eqs. (jA6IIA10|) . with experimental variables oiqL ~ 80 
and f]/(27r) = 10 MHz, so that AR/(27r) = -14.6 kHz, 
rE/(27r) = 105 kHz, and Wg/(27rL) = 16.7 kHz. Because 
the light pulses E and £'* have a finite bandwidth, we 
chose an integration bandwidth of (27r)160 MHz, and we 
have checked that a larger range does not significantly 
affect the results. Solid black curves in Fig. [TCT a-d) plot 
the first approximations in Eqs. (jAlllETs)) : dashed black 
lines show the corresponding Fe — > 00 expressions in Eqs. 



Let us now compute S{z,t). From Eqs. and jl]), we 
have 



S{z,Lj) 



gVNn 
F 



t{z,ijj) — i — 1 (z,w) 

Ahf 



4I6) 



where £{z,uj) and £'*{z,ui) are given in Eq. ()A3[) 
We can then write 



S{z,t)^ J dt'£{0,t^t')hi{z,t') 

+ J dt'£{0,t-t')h2{z,t') 
+ [ dt'£'*i0,t-t')h3{z,t'), 



(A17) 



where hi describes pure EIT, while h2 and are the 
results of FWM. Functions hj can be computed as 
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FIG. 10. (Color online) Graphs of (a) /i(L,t), (a') -hi{L,t), (b) f2{L,t), (b') -h2{L,t), (c) lm[fi{L,t)] = -lm[g3{L,t)\, (c') 
—\m\hz{L,t)\, and (d) g2{L,t). For the calculations, f2/(27r) = 10 MHz, aoL = 80. Red curves show the result of numerical 
integration of the respective expression in Eqs. (jA6HA10l IA18tfA20[l . Solid black curves show the approximate forms of the 
integrals given in Eqs. (jAlHAlSI IA23|) . without taking the limit Fb — ^ oo. The dashed black curves incorporate the Fe — >■ oo 
approximations in Eqs. HA11IIA15I IA23|) . 



h,{z,t') + h,{z,t') = ^ I dc^^^e-^ 



cosh(^z) 

r — 



a 2AR,(r - lu) 



^hf 



cosh(^z) 



Expanding /12 to ©(l/A^^j) and /13 to ©(l/Ahf), the above expressions simplify to 



1 f gVNn (n^e^""^ - F) . 
^ ' ^ 271 J AMiuJ + iTQ)F 



-j^^ ) --h(e^) 



(A18) 
(A19) 
£-"^''. (A20) 

(A21) 
(A22) 



Taking 5 = S,, 70 = 0, 2ia ~ - xi^, and F « 
we have 

h,{z,t')^-^f,{z,t') (A23) 

for j — 1,2,3, where the expressions for fj (z, t') are given 
in Eqs. (|AlHIA13p . Plugging Eq. (|M3|) into Eq. (|CTf)) 
yields an expression that is proportional to the signal field 
£{z,t) in Eq. (jA4p . Thus, under these approximations, 
we obtain Eq. (fT6)) . which is, remarkably, the usual EIT 
relation. Specifically, in the limit of an infinitely wide 



EIT window, S{z, t) can be found by plugging Eq. (|A23|) 
[with the corresponding Fe — ^ 00 expressions for fj from 
Eqs. (|A1HIA13|) ] into Eq. (|JT7)) to yield Eq. (HI]). The 
expressions for \hi{L,t)\, \h2{L,t)\, and |Im[/i3(L, t)] | are 
plotted in Fig. [10] (a', b', and c'), respectively. As before, 
red traces depict the results of numerical integration of 
Eqs. (|A18I - R20)) with the same input parameters as in 
the fj and gj analysis. Solid black curves plot Eq. (|A23[) 
using finite Fe expressions for fj from Eqs. (jAlHIAlS)) . 
Dashed black lines show the corresponding Fe 00 ex- 
pressions. 
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